Competition of multiband superconducting and magnetic order in 
ErNi2B2C observed by Andreev reflection. 
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PACS 74.45.+C - Proximity effects; Andreev effect; SN and SNS junctions 
PACS 74.50.+r - Tunneiing ptienomena; point contacts, weak finks, Josephson effects 
PACS 74 . 70 . Dd - Ternary, quaternary and muftinary compounds (inciuding Ctievref ptiases, boro- 
carbides, etc.) 

Abstract. - Point contacts (PC) Andreev reffection dV/dl spectra for tfie antiferromagnetic 
(Tn — 6K) superconductor (Tc ~11K) ErNi2B2C have been measured for ttie two main crystai- 
iograptiic directions. Observed retention of the Andreev reflection minima in dV/dl up to Tc 
directiy points to unusual superconducting order parameter (OP) vanishing at Tc. Temperature 
dependence of OP was obtained from dV/dl using recent theory of Andreev reflection inciuding 
pair-breaking effect. For the first time existence of a two superconducting OPs in ErNi2B2C is 
shown. A distinct decrease of both OPs as temperature is lowered below Tn is observed. 



Introduction. — The family of quaternary nickel 
borocarbides superconductors i?Ni2B2C, where i? is a 
rare-earth element or Y, has attracted worldwide atten- 
tion both because of a relatively high critical temperature 
Tc, up to 16 K for R—Lu, and especially from the point of 
view of competition between superconducting and mag- 
netic ordered states in the case of i?=Tm, Er, Ho, Dy, 
where energy scales for the antiferromagnetic and super- 
conducting order can be varied over a wide range (see, e. g., 
Refs. [1,2] and further Refs. therein). The compound with 
i?=Er and Tc ~11K is interesting for two reasons [1]: be- 
low (Tn ~6 K) incommensurate antiferromagnetic order 
with spin density wave occurs and weak ferromagnetism 
develops below T^fm — 2K [3]. Both phenomena are, in 
general, antagonistic to superconductivity, so that com- 
petition between superconducting and the magnetic state 
should take place in this compound. Additionally, the su- 
perconducting ground state in borocarbide superconduc- 
tors is expected to have a multiband nature [4, 5] with a 
complex Fermi surface and different contributions to the 
superconducting state by different Fermi surface sheets. 
Therefore, determining the influence of these magnetic 
states on a possible multiband superconducting ground 
state or multiband order parameter (OP) in ErNi2B2C is 



a challenge. 

Previous tunnehng (STM/STS) and point contact (PC) 
spectroscopy results have left some open questions regard- 
ing the coexistence of superconductivity and magnetism 
in ErNi2B2C. STM/STS measurements of [6] show a small 
feature, namely, decreasing of the superconducting gap be- 
low Tn nearly within error bars, which was not reproduced 
in subsequent experiments [7]. Early PC data on polycrys- 
talline samples [9] indicated that the superconducting gap 
has roughly a BCS dependence with only a shallow dip 
around Tn ~6K. Very recent laser-photoemission spec- 
troscopy data [8] show the SC gap decrease (with remark- 
ably large error bars) below the Neel temperature, but, 
at present the laser-photoemission spectroscopy has not 
enough resolution to go deeper in details. 

In this paper we report our detailed directional PC An- 
dreev reflection measurements on single crystal ErNi2B2C 
along the c-axis and in the ab-plane. Our results show 
for the first time the presence of two dominating OPs in 
ErNi2B2C, which differ by a factor of about two, and ap- 
preciable depression of both OPs by the antiferromagnetic 
transition is found. 
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Experimental details. — We have used single crys- 
tals of ErNi2B2C grown by the Ames Laboratory Ni2B 
high-temperature flux growth method [10]. PCs were es- 
tablished both along the c axis and in perpendicular di- 
rection by standard "needle-anvil" methods [11]. The 
ErNi2B2C surface was prepared by chemical etching or 
cleavage as described in [12]. As a counter electrode, edged 
thin Ag wires (0=0. 15 mm) were used to improve mechan- 
ical stability of PCs in comparison to use of a bulk Ag 
piece. We have measured the temperature dependence of 
dV/dI{V) characteristics of such N-S PCs (here N denotes 
a normal metal and S is the superconductor under study) 
in the range between 1.45 K and Tc for several contacts 
oriented both along the c-axis and in the ab-plane. In the 
paper we demonstrate results of analysis of 60 dV/dI{V) 
along the ab-plane measured for the same PC at different 
temperatures between 1.45Kand llKand of 46 dV/dI{V) 
along the c-direction for another PC ^ in the same tem- 
perature range. 

Results and discussion. — To determine the OP 
from the measured differential resistance curves dV/dI{V) 
we used recent theory [13] of Andreev reflection in PC, 
which includes the pair-breaking effect by magnetic impu- 
rities. The last assumption is reasonable, because of the 
presence of the local magnetic moments of Er ions. The 
fit of the measured curves using equations as (1) in [14] 
has been performed. As fit parameters the superconduct- 
ing OP A ^, the pair-breaking parameter 7 = 1/(tsA) 
(here Ts is the spin-fiip scattering time) and the dimen- 
sionless barrier parameter Z have been used. Although 
the dV/dl curves shown in Fig. 1 exhibit one minimum 
for each polarity, as in the case of ordinary one gap su- 
perconductors [11, 15], to fit dV/dl in full we had to use a 
two-OP(gap) approach adding the corresponding con- 
ductivities as made in [12] for LuNi2B2C. The contribu- 
tions of these conductivities account for the part of the 
Fermi surface containing a particular OP. Thus, for the 
two-OP model, experimental curves are fitted ^ by the 



^The PC resistance is 36 Q along the c-axis and 10.5 f! in the 
ab-plane. The PC diameter estimated by Wexler formula (see 
[11], pages 9, 31) is about 7nm and 14nm, respectively, using 
pi ~ 10~^^n cm^ [4]. At the same time a mean free path I is 28 nm, 
using p ~ 3.5~^f2 cm [10] just above Tc. Therefore mentioned PCs 
are close to the ballistic limit d < I. 

■^Assuming that pair breaking is by magnetic impurities, the en- 
ergy gap Aq and the OP A are related as follows: Aq = A(l — 
72/3)3/2 [13], 

^Not only does the one OP approach give a worse fit (especially 
at minima position and at maximum, see insets in Fig. 2) of the 
experimental data such that the rms deviation is 2-3 times higher 
compared to the two OP fit, it also requires a varying Z parame- 
ter. On the contrary, at two OP fit Z parameters remain constant, 
equal to 0.77 and 0.6 for the ab-plane and c-direction, respectively. 
This is important because there is no physical reason for the barrier 
parameter Z to be temperature dependent. 

■^Before the fit, dV/dl curves were normalized to the dV/dl curve 
measured above Tc and symmetrized. The fit was done between 
± 8mV, to avoid contribution from the phonons seen, e.g., as an 
inflection point around 10 mV for some curves in Fig. 1. The fit was 
done in two stages. At first we keep both 71,2 coefficients equal to 
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Fig. 1: (Color online) Raw dV/dl curves of ErNi2B2C-Ag PCs 
in the ab-plane and in the c-direction at indicated tempera- 
tures. For clarity, only several representative curves from the 
total 60 for the ab-plane and 46 for the c-direction measured 
at different temperatures are shown. 



following expression: 



dV 
U 



S 

If (Al, 71,^1) + (1 - X)|f (A2,72,^2) 



(1) 



Here, the coefficient K refiects the contribution of the part 
of the Fermi surface having the OP Ai, S* is the scaling 
factor to match the amplitude of the calculated and the 
experimental curves. 

Before discussing the fitting results, we point out the 
unusual specific behavior of the measured dV/dl ^ (see 



zero. As a result of such fit Ai 2 values shown in Fig. 3 were obtained, 
while variation of Z and K was within 10%. If we will hold K on 
this stage strictly constant, it results only in more scatter (noise) 
for Ai_2, but their overall behavior remains the same. To improve 
the fit on the second stage, we used obtained Ai 2 and variate 71 2 
and K. As a result, obtained K and 71,2 are shown in Figs. 4 and 5. 
After this the theoretical curves were almost indistinguishable from 
the experimental ones (see insets in Fig. 2). 

^Presented in the paper PCs have survived about 36 hours (c-dir) 
and 50 hours (ab-plane) of measurements. The dV/dl temperature 
series for these PCs are the most full, therefore they presented in 
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Fig. 2: (Color online) Reduced position of the minima in 
the raw dV/dl curves for ErNi2B2C and that for LuNi2B2C 
from [12]. Insets: comparison of two OPs (solid line) and one- 
OP (dashed line) fitting of the reduced experimental dV/dl 
(symbols) at 1.45 K. 



Fig. 1). First, the distance between rfV^/dJ-minima shown 
in Fig. 2, which is often taken as a rough estimation of 
the superconducting gap value, increases with tempera- 
ture before decreasing on approaching - quite different 
behavior from the nonmagnetic LuNi2B2C. Second, the 
dl^/d/- minima for ErNi2B2C persist up to temperatures 
close to Tc (see Fig 2, upper panel) though with a small 
amplitude, leading to a supposition of the presence of the 
second OP. From these direct observations, a nontrivial 
behavior of the superconducting OP parameters are ex- 
pected in ErNi2B2C. 

Indeed, from the two-band model fitting both OPs Ai 
(large) and A2 (small) diminish on entering the antiferro- 
magnetic state around 6K (see Fig. 3). Qualitatively the 
same behavior has OP determined by one OP fit (Fig. 3, 
triangles). This is qualitatively consistent with the tem- 
perature dependence of the superconducting gap deter- 
mined by tunneling in [6], by laser-photoemission spec- 



the paper. Of course, there were other PCs with dV/dl of lower 
quality or which did not survive temperature sweep in the whole 
range between 1.45 K and Tc. Nevertheless, there were a few of PCs 
which had dV/dl similar to the presented in the paper, supporting 
our observations. 
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Fig. 3: (Color online) Temperature dependence of the large OP 
Ai (circles), small OP A2 (squares) and OP A determined by 
one OP fit (triangles) for ErNi2B2C for the two main crystal- 
lographic directions. In the upper panel closed (open) symbols 
show OPs determined during increasing (decreasing) tempera- 
ture. The same meaning have closed (open) symbols in Figs. 4 
and 5. 



troscopy data [8] and also with the upper critical field 
[2, 16] and the superconducting coherence length behav- 
ior [17] in the vicinity of T^r. The theories of coexistence 
of superconductivity and antiferromagnetic state also pre- 
dict such OP suppression below Tjv, ® e. g., by antiferro- 
magnetic molecular field [18]. 

Further, the large OP Ai(T) may be described by a 
BCS dependence above Tjv ~ 6K in the paramagnetic 
state with extrapolated T* ~14.5K, close to that of non- 
magnetic i?Ni2B2C (i?=Lu, Y). On the other hand reten- 
tion of the Andreev reflection minima in dV j dl up to Tc 
(see Fig. 2) result in unconventional abrupt Ai(T) van- 
ishing near Tc. Note, that to fit experimental curves, 
not only OP but also large OP relative contribution K 
(see Eq. (1)), must be temperature dependent (see Fig. 4). 
Shown in Fig. 4, the decrease of K with temperature points 



^Various theories of antiferromagnetic superconductors including 
the affect of spin fluctuations, molecular field, and impurities on the 
A behavior (see e. g. H. Chi and A. D. S. Nagi, J. Low Temp. Phys. 
86, 139 (1992) and Refs. therein) in support of our observation will 
be discussed in the forthcoming extended publication. 
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Fig. 4: (Color online) Temperature dependence of the contri- 
bution K (see Eq. 1) of the larger OP to the dV/dl spectra. 
Curves represent a polynomial fit simply to guide the eye. 



to diminution of the " superconducting" part of the Fermi 
surface with the large OP by approaching Tc, which re- 
sults in collapse of the large OP at this point. 

The mentioned decrease of K correlates with behavior 
of the pair-breaking parameter 7 shown in Fig. 5. It ap- 
pears that 71 is always larger than 72 above 2 K, that is the 
pair-breaking effect is stronger in the band with the larger 
OP. This is in line with the conclusion, that the different 
bands are differently affected by magnetic order, made 
in [5, 19] by band structure analysis of the coexistence of 
superconductivity and magnetism in the related antifer- 
romagnetic superconductor DyNi2B2C, i.e., some bands 
provide a basis for superconductivity while others are im- 
portant for the magnetic interactions. Here we should 
add, that the S parameter in (1) has a maximal value of 
0.5 at low temperature, the same value as the normalized 
zero-bias tunnel conductivity obtained by STS study of 
ErNi2B2C [7]. So, both observations suggest that nearly 
half of the total Fermi surface (or bands) is nonsuper- 
conducting in ErNi2B2C. Formation of a superzone gap 
at antiferromagnetic transition seen in transport measure- 
ments [16] may be responsible for this. 

From Fig. 5 it is also seen that 7 has a maxima at tem- 
perature close to Tn and also close to the appearance of 
weak ferromagnetism around 2 K, which is reasonable. At 
both of these temperatures an increase in pair breaking is 
expected due to increasing spin fluctuations accompanying 
the corresponding transitions. 

Conclusion. — This study demonstrates that the 
two-band approximation with two OPs, including pair- 
breaking effects, suits better for describing the PC An- 
dreev reflection spectra in ErNi2B2C pointing for the first 
time to presence of multiband superconducting OP in this 
compound. The values and the temperature dependencies 
of the large and the small OPs have been estimated for the 
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Fig. 5: (Color online) Temperature dependence of the pair- 
breaking parameter 7. Curves represent polynomial fit simply 
to guide the eye. 



ab-plane and in the c- direction. It is found that in the 
paramagnetic state both OPs can be described by BCS 
dependence, but formation of the antiferromagnetic state 
below Tjv — 6K leads to a decrease of both OPs. The 
pair-breaking effect is found to be different for the large 
and the small OP indicating that the different bands are 
affected differently by magnetic order. This may be the 
reason for the observed abrupt vanishing of the larger OP 
at Tc. It is interesting that extrapolation of the larger 
OP by "conventional" BCS behavior above T^r results in 
T* ~14.5K, similar to nonmagnetic YNi2B2C, so that 
2A^°"(0)/kBT;: ~ 4.25 and 4.7 for the ab-plane and c- 
direction, respectively. BCS extrapolation gives for the 
small OP 2A=^<=^(0)/kBTc ~ 4.1 (ab-plane) and 3.5 (c- 
dir), while for the one OP fit 2A="^^(0)/kBTc ^ 4.6 (ab- 
plane) and 4.5 (c-dir) pointing, in general, to moderately 
anisotropic and strongly coupled superconducting state in 
ErNi2B2C. 
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